The Interacting Impurity Josephson Junction: Variational Wavefunctions and Slave 

Boson Mean Field Theory 
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We investigate the Josephson coupling between two superconductors mediated through an infinite 
U Anderson impurity, adapting a variational wavefunction approach which has proved successful 
for the Kondo model. Unlike the Kondo problem, however, a crossing of singlet and doublet state 
energies may be produced by varying the ratio of Kondo energy to superconducting gap, in agreement 
with recent work of Clerk and Ambegaokar. We construct the phase diagram for the junction and 
discuss properties of different phases. In addition, we find the singlet and doublet state energies 
within a slave boson mean field approach. We find the slave boson mean field treatment is unable 
to account for the level crossing. 

PACS numbers; 74.50. +r, 73.40. Gk 



I. INTRODUCTION 

The Ambegaokar-Baratoff formula, /c = 7rA/2ei?, re- 
lates the critical current in a Josephson junction to the 
superconducting gap A and the normal state junction 
resistance R. This result is perturbative in the elec- 
tron tunneling amplitude t; the normal state conduc- 
tance is proportional to |tp when t is small. In certain 
instances, however, it may be that the tunneling is me- 
diated through a magnetic impurity, rather than taking 
place directly from one superconductor to the other. This 
situation has been considered by a number of authors 
10-^. The principal result is that magnetic impurity- 
mediated tunneling results in a negative contribution to 
Ic ■ As Kulik originally argued ^ , the magnetic impurity 
gives rise to an effective spin-flip hopping amplitude tgf 
between the superconductors. Since spin-flip tunneling 
results in a sign change of the Cooper pair singlet (tJ, 
to it) [H, the critical current is Ic oc \t\'^ — |tsfP- When 
Ic < 0, one has a 7r-junction, for which the ground state 
energy is minimized when the phase difference between 
the superconductors is S = tt. Such 7r-junctions break 
time-reversal symmetry (T), and a ring containing a sin- 
gle TT-junction will enclose trapped flux 

This analysis suffices in the perturbative limit where 
t is small. If the impurity level energy is Eq and the 
Coulomb integral is U, the condition for a magnetic 
ground state (and a 7r-junction) is U > — £g > In 
the nonperturbative regime, a new energy scale arises: 
the bare impurity level width T Trpjtp, where t is 
the electrode-impurity hopping matrix element and p the 
electrode density of states. The T = phase diagram as 
a function of — Eq/A, U/A, and F/A was investigated by 
the authors within the Hartree-Fock (HF) approxima- 
tion, where it was found that 7r-junction behavior occurs 
for U > —Eq > 0, provided F is sufflciently small (F < C/). 
However, when A = 0, it is known that the HF approx- 
imation is unable to describe the formation of a Kondo 



singlet at energy scales below = VFexp(— 7r|eo|/2r), 
where W is the half-bandwidth in the electrodes. In our 
problem, then, one expects that the Kondo effect will 
be mitigated whenever A>T^. Roughly speaking, if 
A > Tk, the ground state of the system is a Kramers 
doublet, and breaks time reversal symmetry, whereas if 
Tj^ > A, the ground state is a hybrid singlet formed 
from electrons on the impurity and in the superconduct- 
ing electrodes [|j . 

Recently, Clerk and Ambegaokar |Q applied a general- 
ization of the non-crossing approximation (NCA), a par- 
tial summation scheme, to attack this problem in the 
U ^ oo limit 1^. Adaptation of this method to the in- 
teracting impurity Josephson junction allows one to see a 
transition from 0-junction to 7r-junction when the super- 
conducting gap becomes comparable to the Kondo tem- 
perature. 

In this paper, we also explore the C/ — s- oo limit, using 
two approaches. The first is a variational wave function 
calculation, similar to that used by Varma and Yafet in 
the Kondo problem We generalize this wave func- 

tion in two respects. Firstly, there are two superconduc- 
tors connected to the impurity. This is in contrast with 
the usual setup of Kondo problem where we have only 
metallic electrode. Secondly, the spin-^ state must be 
considered as well. We find, in agreement with ref. 0, 
that a first order transition occurs at T^/A « 1. We 
also show how this transition may be precipitated by a 
change in the phase difference S, as we found in ref. 

The second calculation we describe is the slave boson 
mean field theory. While the this method does confirm 
that the singlet state gives a 0-junction and the Kramers 
doublet a 7r-junction, within the static slave boson mean 
field solution the singlet is always lower in energy (or the 
states are degenerate). The NCA, which goes beyond the 
mean field level, is able to describe the transition. 
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II. VARIATIONAL WAVEFUNCTION 
APPROACH 

Wc start with the grand canonical Hamiltonian K 

H-iiN, 



a q ^ 



(1) 



where a labels the electrode. S^qa is the dispersion in the 
a electrode relative to the chemical potential, is the 
modulus of the superconducting gap and Sa its phase, 
Afa is the number of unit cells in electrode a, ta is the 
hopping amplitude from electrode a to the impurity, and 
£q is the bare impurity energy. We set U — » oo, which 
leads to the constraint "^^cl-c^ < 1. 

Defining the angle 6qa = ta.n^^{Aa/£,qa) and the 
usual BCS coherence factors Wq„ = cos{^6qa), Vqa = 
sm{^9qa) cxp(iSa), we express /C = /Cq + /Ci in terms of 
the Bogoliubov quasiparticle operators: 



q,a 



q,a ^ 

+^9a(7qaT4 " 7qai4) + ^Qa(7qaiC^ " Iqat^'O 



where i;^^ = y^Cqa + A^, and = l^al/V^- 

We now two variational many-body states for the U = 
oo limit: a singlet, 



S) 



t .t 



■ Tqat4) 



1 ' 

9'.' 



(2) 



and a doublet, 



(3) 



where | ) is the fermion vacuum (the other doublet state 
I D J, ) is obtained by rotating the spins by tt about the 
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FIG. 1. Singlet (solid) and doublet (dot-dash) energies 
versus phase difference for Eq — —1.6 (upper left), e^-, — —1.8 
(upper right), Sq = —1.9 (lower left), and Sq = —2.0 (lower 
right). In all cases A = 1 and = Fj^ = O.Stt. 



2/-axis). Here, C^',^ = C^^', , C^,^ = C^"/, and b<^,^ 
—Dqq, . We next set to zero the variations 



5{<i\Ko + Ki\^) - E5{-^\-^) = Q 



(4) 



where | \1/ ) is | S ) or | D t ) , to obtain equations relating 
the variational coefficients. We find that A and the ma- 
trix Cqq, may be expressed in terms of the coefficients 

Bqa, and similarly A and the matrices Cqq! and Dqq, 

may be expressed in terms of the coefficients Bq^- We 
then obtain the two eigenvalue equations for the singlet 
and doublet energies E and E, respectively: 



E 

q'a' 



and 



E 

q'a' 



2 Vqa Vq'a' tqa tq,^, 



^qa '^q'a' ^9" iq' a' 



E 



^ ~ ^qa ~ ^> 



E 



0-E 



^qa 

ul,^, t 



Eq'a' 



Bq'a' 



q'a' ^q'a' 



Uqa '^q'a' ^Qa ^q'a' 



q',a' ^ ^qa ^q'a' 



Vqa Vqf^t tqa tq' a' 



B. 



qoi 



(5) 



E- 



E-e^ 



E. 



qa 



^q'a' 



Bq'a' 



E-E ^ '^\Vq'a'?t\>a' 



q',a' ^ ^0 ^qa ^q'. 



B. 



qa 



(6) 



We solve these equations numerically for the symmet- 
ric case Al = Aj^ = A, = = t, Tj^=T^ = F. The 
normal state of each electrode is described by a flat band 
of width 2W; we use W/A = 10 in our calculations, but 
the general features are rather insensitive to the value of 
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A = 2Wexp(-7r|eo|/2r) = 2T^ 



(9) 



FIG. 2. Variational wavefunction phase diagram for the 
symmetric case = Fp^ = F. 



W. Energy versus phase difFerence, E{5), is plotted in 
figure |l| for A = 1, F = O.Stt, for four different values of 
— Eq. When Eq = —1.6, the singlet state is the ground 
state for all (5, while for Eq = —2.0, the ground state is 
always the doublet. These are 0- and 7r-junctions, respec- 
tively. For Eq = —1.8, the curves cross, and the ground 
state energy has a kink as a function of 6. Both (5 = and 
(5 = TT are local minima in E{5), with S — the global 
minimum. Using the terminology of ref. , this is a 0'- 
junction. The final case of the 7r'-junction is reflected in 
the curves for Eq = —1.9, where ^ = is a local minimum 
and 6 = 7: the global minimum of the energy. 

The phase diagram is displayed in fig. 2. In the HF 
treatment ||^, one also finds a transition from tt to tt' 
to 0' to 0-junction phases as F is increased from zero. 
However, the critical value of F there turns out to be 
inversely proportional to In U. There are two reasons 
for this: first, the HF calculations assumed an infinite 
bandwidth, and second, HF is unable to describe the for- 
mation of a Kondo singlet. Here, the transition from 
0-junction to 7r-junction occurs roughly for A ^ Tj^, in 
agreement with ref. 0. Indeed, when < — Eq,F <C A 
this relation can be derived within perturbation theory. 
The perturbative value of energy of the singlet state is 



En 



4F, 2W 
In-r-, 

TT A 



while the energy of the doublet is 
E 



2F, 2W 
— m — — . 

TT A 



(7) 



(8) 



Comparing these two we find that the symmetry of the 
ground state changes when 



This predicts a straight line for F versus — Eq with a slope 
F/(-Eo) = 7r/21n(2VF/A). With W = lOA, the slope is 
0.524, and the line would be bounded by the — 0' and 
TT — tt' phase boundaries in fig. 2 (it is not shown for the 
purposes of clarity). 



III. SLAVE BOSON MEAN FIELD THEORY 

We consider an extension of the model of eqn. (|l|) by 
introducing a flavor -1 m which runs from 1 to A^f. The 
Hamiltonian is then 



/C = 



■00 



r + 



EEE 

a m—1 q 



E (^9" + /^bS^) 



E E 

m=l a=1,l 



(Eq + Mb^ct). 



'maimer J 



(10) 



where is the Zeeman energy. The impurity level 

occupancy satisfies the constraint J^m a ^^Ina'^ma — ^■ 
We refer to this as model I. In a slightly different 
large- iVf extension (model II), we rescale the hopping 
amplitudes —* ta/V^^f and write the constraint as 

Em,<7 4naCrr,a < ^Nf. In both CaSCS, < T < 2. 

Introducing a slave boson b and a Lagrange multiplier 
A to impose the constraint, we evaluate the impurity con- 
tribution to the free energy at the mean field level, as- 
suming both b and A to be static. The impurity free 
energy per fiavor is then 

Fimp/Nf = e-fiBB+ p{e - Eo)(|6|' - r) 

OO 

f{u; + HsB) Im In H{uj + iO+) 



2|6|2F„Asin(i(5) 



A2 



VA2 



(11) 



where F^^ = Fg = \/rLrR ^'^^"^ = T:pc\ta^ [pa 
is the bare density of states per unit cell in electrode 
a), 6 = — is the phase difference between the two 
superconducting electrodes (we assume Ai = A2 = A), 
and f{uj) = [exp(w/T) -I- 1]"'^ is the Fermi function. The 
Lagrange multiplier A is absorbed into the renormalized 



impurity level energy, e 



Eq + A. 



The factor p in the 
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second term is 1/Nf in model I, while p = 1 in model 
II. Hence it is model II which generates a true large- TVf 
expansion, u with all terms in the impurity free energy 
^imp proportional to TVf . 

From H{Q) < 0, H{A-) = +00, and dH/dio > 0, we 
conclude that there is a unique solution to the equation 
H{u!) = on the interval u; e [0, A]. Call this root f2. 
We use the external field i? as a Lagrange multiplier so 
that we may fix the total value of S^. Making a Legendre 
transformation to G-^^^iS'') = F.^^{B) - 2/^B-BS'^ with 

= — i,0,-|-i, we set dF-^^^/d{^^B) — 0, and obtain, 
at T = (assuming < ^b-B < A), 



G?„,p/7Vf-e + p(£-£o)(|6p-r)+^ 
1 



(12) 



A^- dcjImlniJ(cj-iO+)-r2(5s-o , (13) 

TI" J A 



where oJq — VW^ + A^; W is the half-bandwidth in the 
electrodes. The mean field equations, obtained by setting 
5GP„p/9|6|2 = and dGl^^^/de = are 



-Pi\b\' 



q) 



dA 
dA 



and the Josephson current is 



2edA 
h dS 



(14) 
(15) 

(16) 



Thus, for < /ie-B < 17, the ground state has — 0, 
while for fl < l/iB^I < A, the ground state has 5^ = 
-isgn(B). 

For p = r = 1, we can show that the ground state 
energies satisfy AGP,„„ = GP„„ — GP„„ > 0, which means 

*— ' lllip lllip lllip ' 

that the mean field slave boson theory cannot describe 
the TT or tt' phases. The energy difference AG?,-jp is min- 
imized at 5 = TT. The value of AG?j^p is an increasing 

function of the bare impurity level energy Sq. However, 
there are no solutions to the mean field equations for 
Eg < e™n ^ _47r-ir^cosh"\wo/A) - 2TJs'',o- Rather, 
the endpoint solution e = |6p = holds. Thus, the best 
we can do is AG^j^^p = 0, but in this case G^^^ = Eq, 
independent of S, and there is no Josephson current. 
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FIG. 3. Mean field parameters versus Sq. A = 1, 
= 0.650, Fg = 0.630, B = 25, 5 = 0, and p = r = 1. 
Solid curve is 5^ = 0, dot-dash is S'^ = ±|. 
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FIG. 4. Mean field parameters and Josephson current 
versus 5 = for £(, = -2, A = 1, F^ = 0.650, Fg = 0.630, 
B = 25, 5 = 0, p = r = 1. Solid curve is = 0, dot-dash is 
5^ = ±i. 



IV. LARGE A LIMIT 



The case A ^ 00 (with U finite) may be solved exactly. 
We begin with the Hamiltonian of eqn. (|^), integrating 
out the fermion degrees of freedom in the superconduc- 
tors [o|. This generates an induced action. 



5i 



ind 



E 




EJ- a I iljJm Ae 



Juj'^ +A2 



(17) 



If the dynamics all occur on frequency scales w ^ A, we 
may ignore the Matsubara frequencies LOm in comparison 
with A. Adding the induced action to the bare action 
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for the impurity, we find the resultant action is that for 
a Hamiltonian 

(18) 

where 

(19) 



VI. BIBLIOGRAPHY 



The ground state energy is Ej^ = Eg for the Kramers 
doublet, and 



for the singlet. Thus, for U < Uc, where 



c/c = 4Wr2-r2sin^(-5) , 



(20) 



(21) 



the Coulomb repulsion is too weak to overcome hy- 
bridization effects, and the ground state is a singlet for 
all Eg. For U > Uc, the ground state will be a doublet 
provided 



-U- _ <e„<-U+ 



(22) 



This allows the possibility of a level crossing as a function 
of S ii < < E^. However, the doublet energy is in- 
dependent of 6 in this model, hence there is no Josephson 
coupling in the doublet ground state. 
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V. CONCLUSIONS 



Using a generalization of the variational wavefunc- 
tions applied in the study of the Kondo effect |^,^, we 
have demonstrated that the U = oo interacting impurity 
Josephson junction exhibits a first order phase transition 
for A « Tj^. When A<T^, the ground state of the sys- 
tem is a singlet, and E(5) is minimized at S — and 
maximized a,t S — n. When A>T'j^, the ground state is 
a Kramers doublet, the stability of (5 = and ^ = tt is re- 
versed, and the system forms a 7r-j unction. For A ~ T^, 
the singlet and doublet energy curves may cross, in which 
case the ground state energy has a kink as a function of 
5 ^ and is strongly nonsinusoidal. 

We also solved for the junction's properties within a 
slave boson mean field theory. Unfortunately, this ap- 
proach is unable to identify a phase transition, and the 
singlet state always is lower in energy than the doublet. 
One must go beyond the static mean field approach, as 
has recently been accomplished in ref . Q , to see the tran- 
sition. 
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